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Abstract 

We extend a recent theory of parametric correlations in the spectrum of random 
matrices to study the response to an external perturbation of eigenvalues near the 
soft edge of the support. We demonstrate by explicit non-perturbative calculation 
that, for systems with unitary symmetry, the two-point function for level density 
fluctuations becomes, after appropriate rescaling, a universal expression. 
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Random matrix theory (RMT) [1] has recently been the focus of much attention 
because of its wide range of applications in many branches of physics, such as com- 
plex nuclei [2], disordered metallic grains [3], quantum chaotic systems [4], random 
surfaces [5], disordered conductors [6], QCD [7] and more recently strongly corre- 
lated many-particle systems [8] . One of the most striking and robust predictions of 
RMT is the celebrated Wigner-Dyson statistics [9] , which states that the probability 
distribution of level spacings is a universal function if the eigenvalues are measured 
in units of the mean level spacing. Another important feature of RMT is that all 
physical systems are classified according to its symmetries in just three universality 
classes, namely: unitary, orthogonal and symplectic. 

A significant extension of RMT has recently been put forward by Simons, Alt- 
shuler, Lee and Szafer [10-11]. In particular, they have shown that in order to 
describe the statistical properties of the dispersion of the energy levels of a particle 
diffusing through a disordered metallic grain in response to an external adiabatic 
perturbation, one has to introduce a new parameter in the theory in addition to the 
mean level spacing. This new parameter has been given the meaning of a generalized 
conductance because of its relation to the rate of energy dissipation of the system 
due to the presence of the external perturbation. The most remarkable result of 
their work, however, is that the n-point correlation function for density of states 
fluctuations becomes a system independent universal expression if the energy levels 
are measured in units of the mean level spacing and if the perturbation parameter is 
renormalized by the mean square gradient of the levels. This striking universality has 
been analytically derived for disordered metallic systems with orthogonal, unitary 
and symplectic symmetries and has been verified numerically in a number of chaotic 
[10,12] and also interacting systems [13]. More recently, Simons, Lee and Altshuler 
[14] have demonstrated that the non-linear a-model provides a common mathemat- 
ical structure in which a number of problems in different branches of theoretical 
physics can be interrelated. 

The universality obtained by Simons et al, and similarly that of Wigner-Dyson 
statistics, applies only for levels away from the edge of the support of the spectrum. 
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This statement can, in principle, be formally justified by means of renormalization 
group arguments. A simple renormalization group procedure has recently been de- 
vised by Brezin and Zinn- Justin [15]. These authors have shown that there are two 
kinds of fixed points in the renormalization group equations: a stable gaussian one 
governing the behavior of the system at the bulk of the spectrum and an unstable 
one governing a small crossover region around the endpoint of the support where 
the average density goes to zero as a power law. The attractive gaussian fixed point 
supports the general validity of the Wigner-Dyson statistics everywhere at the bulk 
of the spectrum. It is therefore natural to expect that the results of refs. [10-14] are 
of similar general validity, although an explicit proof is not yet available. 

At the edge of the spectrum, however, Wigner-Dyson statistics breaks down and 
an entirely new regime of parametric correlations emerges. It is well known that the 
average eigenvalue densities of random matrix ensembles exhibit two kinds of edges: 
a hard edge and a soft edge, where the eigenvalue integration range is bounded and 
unbounded respectively. The hard edge of the Laguerre ensemble is important in the 
description of many physical systems such as disordered metallic conductors [6] and 
QCD [7] and has been the subject of a recent paper [16], in which the universality of 
the parameter dependent two-point level density autocorrelator has been explicitly 
demonstrated. The objective of the present work is to study the response to an 
external perturbation of the eigenvalues of a random matrix near the soft edge of 
the support of the spectrum. We demonstrate explicitly that, for systems with 
unitary symmetry, the two-point correlation function for level density fluctuations 
becomes a universal function after rescaling the levels by an appropriate system 
dependent constant and the perturbation parameter by the mean square gradient of 
the levels. In order to emphasize the robustness of the final result we consider two 
different kinds of unitary ensembles: the Gaussian unitary ensemble (GUE) and the 
Laguerre unitary ensemble (LUE). Physically, they describe systems with broken 
time-reversal symmetry. 

The problem of the behavior of random matrix models near the soft edge of the 
spectrum is of wide interest and has been motivated by the application of RMT to 
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the physics of random surfaces [5]. It has been demonstrated [17,18] that the double 
scaling limit of these models, in which the coupling constant of the theory approaches 
a critical value and the size of the random matrix tends to infinity, is characterized by 
a universal crossover function of scaling variables, whose definitions depend on the 
order of multicriticality of the problem. The crossover region has been shown to be 
around the soft edge of the support of the spectrum and its universality is controlled 
by an unstable fixed point in the renormalization group flow [15]. We remark that 
the two random matrix ensembles studied in the present work have the same order 
of multicriticality and therefore must exhibit the same universal behavior, as will be 
demonstrated by explicit calculation. 

An interesting application of the model we consider can be found in ref. [19], 
in which the string theory of matter with central charge equal one coupled with 
2D quantum gravity is represented as a continuous matrix model similar to the one 
studied in ref. [14], where its connection to the non-linear cr-model and the Dyson 
Brownian motion model [20] has been demonstrated. 

In this work, however, we shall, for the sake of generality, consider as our start- 
ing point directly the Dyson Brownian motion model rather than the continuous 
matrix model of refs. [14] and [19]. The Dyson Brownian motion model describes 
the evolution of the probability distribution P M ({xi},r) of N classical particles at 
positions {xi(r)} moving in a viscous fluid with friction coefficient 7 _1 at tempera- 
ture 
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for the Laguerre ensemble [16], where < Xi < oo. The parameters c and a 
are system dependent constants and (3 = 1,2 and 4 for systems with orthogonal, 
unitary and symplectic symmetries. The fictitious time r can have several meanings 
depending on the particular physical problem. It can, for instance, be the one spatial 
dimension of the target space of the continuous matrix model representing the string 
theory of ref. [19] or it can be related to the perturbation parameter 5u of refs. [10- 
14,21] by r = 5u 2 . The average level density p^(x) = S(x — Xi(u))) eq , where 
(■ ■ -)eq stands for an arbitrary average calculated with the equilibrium distribution 
pv({xi}, oo), can be calculated explicitly by direct integration of Eq. (1). We find 

G (x) _ { (2c l /yn)VN^x^ for |x| < (N/c) 1 / 2 

9 {X) ~\0 for |x| > (N/c) 1 / 2 {2a) 

which is the Wigner semicircle law [1] and 

n L M _ f (c^/tt) V4iV - or 2 for < x < 2{N/c) 1 / 2 

P [X) ~\0 for x > 2(iV/c)V2 ' W 

which is a quarter circle as was first shown by Bronk [23]. Note that p^{x) has a soft 
edge at Xq = y/N/c and Xq = 2 a/ N/c. Near this edge, the average level density 
behaves as p^(x) ~ A^ fc in a region of size ~ N~ k , in which k is a multicritical 
exponent [18], which can be shown to be equal 1/6 for both GUE and LUE. The 
double scaling regime of these ensembles is obtained simply by moving the origin of 
the spectrum to Xq and measuring the eigenvalue spacing in units of order iV -1 / 6 . 
An explicit calculation of the eigenvalue density crossover function can be found in 
refs. [18] and [24]. 

In what follows we shall be interested in studying the two-point function for 
level density fluctuations defined as 

S"(x, y, Su) = S ( x ~ Xi(u))S(y - x 3 (u + 6u))) ^ - p fl (x)p^y) (3) 

in the limit x — > Xq and y — > Xq . 

It has been shown [11,16,21,22] that Eq. (1) can be exactly mapped onto a 
Schrodinger equation in imaginary time. The case (3 = 2, for systems with unitary 



symmetry, is particularly simple since the two-body interaction term vanishes and 
we end up with a theory of free fermions. In second quantization the resulting 
Hamiltonian reads 

oo 

T~i ^ £ p c p c pl (4) 

where e p = {2c/^)(p + 1/2), e p = (2c/ r y)(p + (a + l)/2) and c|, (c p ) creates (an- 
nihilates) a fermion with quantum number p. It is convenient at this stage to 
introduce the field operators ip^(x) and ^(x), where ip^(x) = J2 p 0^ (x)c p and 

*I>U X ) = E P <^<»4' in which 

/o \ 1/4 

0G (a .) = / _£j (2 p pl)- 1 / 2 e- cx2 H p (xV2^) (5a) 

and 

= (r5TTTJo) *- M/> «~ v %"(« s >- (») 

The functions H p (x) and (x) are the conventional Hermite and associate Laguerre 
polynomials respectively. Since 7i M is a free electron Hamiltonian we can use Wick's 
theorem to evaluate S^(x,y, Su). We find 

oo JV-1 

y, <5«) = £ ^(*)^(y) £ ^(x)^(y) exp(2co W 2 (g - p)/ 7 ). (6) 

p=0 q=0 

For large A" and x — > we can use the asymptotics [25] 

e~ x2 / 2 H N (x) ~ 7r 1 /4 2 ^/ 2 + 1 /4 (iV !)i/2 i v-i/i2 Ai((;r _ v W + T)iV 1 /62 1 /2) (7a) 

and 

e-* 2 / 2 ^ 2 ) ~ (-1)^2— 1 / 3 iV" 1 / 3 Az((x - 2iV 1 /2) 2 2/3 i vi/6 )7 (76) 

in which is the Airy function. Inserting (7) and (5) into (6) and taking A" — > oo 

and c -> 0, such that 2c 1 / 2 A^ 1 / 6 -> ^ and 2 2 / 3 c 1 / 2 N 1 /6 _^ ^ we get 
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S^X, y, 8u) =pf dz dz'Ai(p%X - z)Ai(p l Y - z)Ai(p l X + z' 
Jo Jo 

Ai{p»Y + z') eM- Ar^i* + (8a) 
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where X = x — Xq and Y = y — y^. The parameter p$ is related to the average level 
density through = pg Jo°° dzAi 2 (pQX + z). Note that if we make the rescalings 

[10-14] X = p%X, Y = p%Y, Su 2 = Su 2 pfh and S(X,Y,Su) = S^(X,Y, Su) / pf 
we find that 

POO POO 

S(X 7 Y,Su) = dz dz'Ai(X - z)Ai(Y - z)Ai(X + z) 
Jo Jo 

Su 2 

Ai(Y + z') exp(-^-(z + z')) (86) 

becomes a function that does not depend on any physical parameter of the system 
and therefore is a universal crossover function characterizing the response to an 
arbitrary external perturbation of the eigenvalues of GUE and LUE near the soft 
edge of the support of the spectrum. This is the central result of the present paper. 

We can easily show that Eq. (8a) reproduces known results for Su — if we use 
the relation [26] dzAi(x + z)Ai(y + z) = 8{x — y) to write S^(x, y, 0) as 

S»(x, y, 0) = 8(x - y)p^(x) - (K»(x, y)f , (9a) 

where 

K"(z, y) = (At(p^x)At'(p^y) - Ai(f%y)A?(f$z))/(x - y). (9b) 

Equation (9b) constitutes the Airy kernel [24], which is known to characterize the 
soft edge of the GUE and LUE. Equation (9a) is also well known in RMT and is 
usually derived using the orthogonal polynomial technique [1] . 

An important application of Eq. (8a) is to calculate the density of level gradients 
defined as F(v) = (S(v — dx/du)). One can show that [27] 

Su 

F(v)= lim — — S^(x,x + vSu,Su). (10) 
Asymptotic evaluation of (8a) for Su <C 1 yields 

Inserting (11) into (10) gives 



This result is known to be valid at the bulk of the GUE and at the bulk and hard 
edge of the LUE. The present result, valid at the soft edge of both the GUE and 
LUE, leads to the striking conclusion that the density of level gradients is gaussian 
everywhere in the spectrum. This allow us to introduce, also at the edge of the 
spectrum, the concept of a generalized conductance characterizing the rate of energy 
dissipation of the system in response to the external perturbation. It is defined as 
[10] 

0) = p{fV> = p(fV7- (13) 

If we now compare (13) with the rescaling of the perturbation parameter Su that 
leaves the two-point function for level density fluctuations universal, we can see that 
the rescaling can be written as Su = C^ 2 5u, which is exactly the same as that used 
in refs. [10-14]. The rescalings of X and Y differ from that of the bulk because there 
is no local translation symmetry at the edge of the spectrum. 

In conclusion, we have studied the response to an external perturbation of the 
eigenvalues of GUE and LUE near the soft edge of the support of the spectrum. The 
main effect of the external perturbation is to cause the levels to disperse in disjoint 
manifolds, as a result of eigenvalue repulsion or non-integrability of the underlying 
dynamics. We have shown that the two-point level density auto correlator becomes 
a universal function when the levels are measured in units of an appropriate system 
dependent constant and the perturbation parameter is rescaled by the mean square 
gradient of the levels. The resulting universal expression is identical for GUE and 
LUE and describes the crossover region where the average level density goes to zero 
as a power law. In addition, we have demonstrated that the density of level gradients 
is gaussian everywhere in the spectrum. 
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